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Abstract 

Morphologies of moving amoebae are categorized into two types. One is the "neu- 
trophil" type in which the long axis of cell roughly coincides with its moving direc- 
tion. This type of cell extends a leading edge at the front and retracts a narrow tail 
at the rear, whose shape has been often drawn as a typical amoeba in textbooks. 
The other one is the "keratocyte" type with widespread lamellipodia along the front 
side arc. Short axis of cell in this type roughly coincides with its moving direction. 
In order to understand what kind of molecular feature causes conversion between 
two types of morphologies, and how two typical morphologies are maintained, a 
mathematical model of amoebic cells is developed. This model describes movement 
of cell and intracellular reactions of activator, inhibitor and actin filaments in a uni- 
fied way. It is found that the producing rate of activator is a key factor of conversion 
between two types. This model also explains the observed data that the keratocye 
type cells tend to rapidly move along a straight line. The neutrophil type cells move 
along a straight line when the moving velocity is small, but they show fluctuated 
motions deviating from a line when they move as fast as the keratocye type cells. 
Efficient energy consumption in the neutrophil type cells is predicted. 
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1 Introduction 



Amoebic cells have been use d as model sy stems to analyze the molecular 
mechanism of cell locomotion (IPollardl . 2003j). Chemical or mechanical signals 
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arriving on amoebic cell membrane induce the intracellular signaling cascade, 
resulting in the activatio n of Arp2/3 complex to in itiate actin polymerization 
at the front edge of cell ( Pollard and Borisv . 20031 ). Thus elongated actin fil- 
ament drives protrusion at the cell anterior. Myosin II seems t o play roles in 



the retraction process at the cell posterior ( Eddv et al. . 2000l ). and the cell 



moves forward when protrusion at the front is followed by retraction at the 
rear. Although those molecular details have begun to be elucidated, there still 
remain large mysteries on how biochemical reactions drive the coordinated 
change of morphology. 

Morphologies of moving amoebae seem to be categorized into two types. One 
is the "neutrophil" type in which the long axis of cell roughly coincides with 
its moving direction. This type of cell extends a leading edge at the front 
and retracts a narrow tail, whose shape has been often drawn as a typical 
amoeba in textbooks. The other one is the "keratocyte" type with widespread 
lamellipodia along the front side arc. Short axis of cell in this type roughly 
coincides with its moving direction. Fish or amphibian epidermal keratocytes 
keep their stable half-moon shape during locomotion and rapidly crawl along 
the relatively straight line. This stability in morphology and motility has al- 
lowed detailed investigations on keratocytes ( Lee et al. ~ ll993l ). Although there 



is further variety in cell shapes depending on the strength and frequency of pro- 
trusion and structures of cytoskeltons, we here focus on the difference whether 
the shape is head-tail long or laterally long and call the former "neutrophil 
type " and the latter "keratocyte type " . 

Interesting observations are that changes in a single or a small number of 
genes can convert cells from the neutrophil type to the keratocyte type: Breast 
adenocarcinoma cells overexpressing a Ca s-family protei n show ed morpho- 



logical conversion to the keratocyte type ([Fashena et al. . 12002ft . A cell line 



of human mic r ovascular endothelial cell moved in a keratocyte-like manner 
( Kiosses et al. . 19991 : Fischer et al. . 2003| ). An interesting example of conver- 
sion is amiB- Null Dictyostelium cells whose amiB gene is disrupted. The wild 
type Dictyostelium cells are representative neutrophil type cells and move to- 
ward the aggregation center with repeated dynamical protrusion and retrac- 
tion when they are starved. The starved amiB~ cells, o n the other hand, m ove 



unidirectionally with maintaining a half-moon shape ([Asano et all 120041 ) . A 



role of AmiB in motility has not yet been known but this example clearly 
showed that the change in a single gene caused distinct conversion of mor- 
phology from the neutrophil type to the keratocyte type. 



Here, in this paper we ask questions of what kind of molecular feature causes 
conversion between two types of morphologies, and how two typical morpholo- 
gies are maintained. In order to answer these questions, theoretical model is 
necessary to be developed, which can treat motion, sensing, and chemical 
reactions in a unified way. Such mathematical model of amoeba has been 
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developed by iNishimura and Sasai (2004, l2005al lbl). and it has been shown 
that the time lag between the shape change and the intracellular molecular 
diffusion induces cell locomotion with "inertia" -like features: The model can 
explain experimental observations that once cells are boost ed in one direction 



cells can keep moving forward without chemoattractants (iVerkhqvskv et al 
1999) or even toward the direction of decreasing chemoattractants ([jeon et al 



2002). It was predicted that in a "centripetal" distribution of chemoattractants 



cells can show a rotatory m otion by keeping fictitious angular momentum 
( Nishimura and Sasai . 2005b| ). These fictitious inertia or angular momentum 
results from the collective dynamics of morphological change and chemical 
reactions and diffusion. In the present paper we use this model to explain the 
dramatic change of morphology between the neutrophil type and the kerato- 
cyte type. 



2 Model 



As in fluid or solid mechanics, the movement of an extended object can be 
modeled with either way of descriptions, the Lagrange description or the Euler 
description. In the Lagrange description a cell is regarded to be composed of 
many small pieces and equations of motions of these pieces are considered. 
What is advantageous in the Lagrange description is the easiness to treat 
forces and the mechanical balance among them in cell. In the Euler description, 
on the other hand, the movement of cell is described from the coordinate 
fixed in space. This description is useful to treat chemical distributions and 
reactions in cell. In models of amoebae both two descriptions were adopted in a 
hybrid manner (|Bottino et all 120021: iRubinstein et alll2005| ) to treat the force 
balance and the chemical reactions at the same time. Another way to express 
amoebae is to solely use the Euler description. In this case phenomenological 
rules are used instead of treating the for ce balance in equations of motions 
f)Satvanaravana and Baumgaertnerl . 120041 1. The advantage of this approach 
is the simplicity in modeling, which allows highly efficient simulation of cell 
movement. We here adopt the Euler description to calculate a large number 
of trajectories for statistical sampling. 



We introduce hexagonal grids that indicate the simulation space. A cellular 
body is thought of as a set of single connected grids in the grid space. We 
call these grids "cellular" grids. We call remaining grids in the model space 
"external" grids. It has been known that the moving cell repeats adhesion and 
de-adhesion to substratum, suggesting that t he cel l movement is regulated by 
a "clutch-like" mechanism ( Smilenov et al. . 1999} ) . Here, we simply assume 
that the cellular grids are linked to the substratum by adhesive molecules and 
do not slide. Although there are complex biochemical pathways to activate 
and inactivate adhesive molecules such as integrins and associated proteins, 
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we omit modeling these processes for simplicity. In our model conversion of 
the external grid to the cellular grid and conversion of the cellular grid to the 
external grid are assumed to represent adhesion and de-adhesion processes, 
respectively. If a cellular grid connects at least one external grid, we call this 
grid a "membrane" grid. If a cellular grid is not a membrane grid, we call 
that grid "cytoplasm" grid. There are three different molecular densities de- 
fined in each cellular grid: activator, A, inhibitor, /, and actin filaments, F. In 
addition, concentration of the external chemical signal accepted by receptors, 
S, is defined at each membrane grid. Activator and inhibitor are regulating 
proteins. Though the regulating biochemical network has not yet been fully 
identified, PI3K and PTEN are candidates of prot e ins w orking as activator 
and inhibitor, respectively (jLevchenko and Ig lesiad . l2002h . It is known that 



PI3K phosphorylates PI(4,5)P2 to PI(3,4,5)P3, whereas PTEN dephosphry 



lates PI(3,4,5)P 3 to PI(4,5)P 2 , working in an antagonistic way (Iiii ma et al 
200|). 



We use "rule based dynamics" in which several rules are called randomly. The 
following paragraphs explain those rules. 

(1) Kinetics: Both a ctivator and inhibitor are pro duced by the stimulation 
of the external signal ([Levchenko and I glesiasl . l2002h . The activator enhances 



polymerization of actins, whereas the inhibitor suppresses the polymerization. 
First, this rule selects a grid in the cellular domain randomly. If densities of 
activator, inhibitor and actin filaments at the selected grid j are expressed as 
Aj, Ij and Fj, respectively, those variables are changed obeying the following 
equations: 



Aj = Aj + aSj - k a Aj (1) 
I'j = I, + (3S 3 - k p I 3 (2) 

[ 7 - kfF, > h) 
F( = Fj + \ f 3 y> , (3) 

I —kfFj (otherwise) 

where a, j3, 7, k a , hp, kf, 7 and h are constants, a and (3 are the creation 
rates of activator and inhibitor, respectively, and 7 is the rate of polymeriza- 
tion of actin filaments. k a and kp are the degradation rates of activator and 
inhibitor, and kf is the de-polymerization rate of actin filaments. Increase in 
h decreases the possibility of polymerization of actin filaments. Sj indicates 
the concentration of chemoattractants or the strength of the external signal 
at the jth grid. Sj is set to zero if the jth grid is in cytoplasm. 



(2) D iffusion: Only the inhibitor diffuses into the whole cytoplasm ( Levchenko and Ielesiasl . 



2002| ) . This rule selects a grid from the whole cellular domain. At the selected 



jth grid and its nearest cellular Zth grid, L and l\ obey the following equations: 
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I' = h - DL 



T 4- ^ 
J; H ■ 



(4) 
(5) 



where D is the diffusion constant, n is the number of the nearest cellular grids. 
D should be smaller than 1 by definition. Since activator and actin filaments 
do not diffuse and the external signal Sj is nonzero only at membrane grids, 
activator and actin filaments distribute around membrane grids. 

(3) Keeping the cell: We also give a rule to prevent cell from breaking 
into pieces. The cellular volume is kept and the cellular surface length is con- 
strained to be as small as possible. This rule randomly selects a grid from 
the membrane. Then the rule decides either to remove the grid or to add a 
new cellular grid around the grid. This rule checks the cellular tension by 
calculating a cost function, E, as E = (V — Vq) 2 + cL 2 , where V is the total 
number of cellular grids and L is the total number of membrane grids, and 
Vq and c are constants. Vq is the equilibrium number of cellular grids and c 
is a " stiffness" -like factor. When E' denotes the cost function after a cellular 
grid is removed or added, we define the probability P as P = exp f~ ^f^ ); 
where kT is a constant to control the extent of fluctuation. We generate a 
random number between and 1 and then compare the number with P. If 
the number is smaller than P, we "undo" the event of removing/adding. Note 
that if "removing" is chosen, the values of A, I and F in the removed grid 
are added into the nearest cellular grid. If "adding" is chosen, the three values 
are inherited from a randomly selected nearest grids, and three values of the 
selected nearest grid are reset to zero. 



This rule is related to contraction of the rear of the cell. Recent experimental 
studies have elucidated that interactions among myosin II and actin filaments 
generate forces to contract the rear. For cellular locomotion, protrusion at the 
leading edge seems to h ave a prime ro le of cellular "decision making" about 
where the cell should go (Weber, 2006t ). and the myosin-generated contraction 
works to weaken the load of protrusion at the leading edge. In this way the 
rate of contraction may be regulated so that the cell size is kept as constant 
as possible. We should note that not only the myosin II contraction but also 
the cellular cortical tension exerted by actin-myosin network is a factor of 
keeping the cell si ze, where myosin I seems to contribute to the generation 



of cortical tension ( Dai et al. . 19991 ). Thus, the effects of the rear contraction 



and the cortial tension could be modeled by imposing the rule of keeping the 
cell size. Several theoretical studies on cell locomotion explicit ly or implicitly 



assumed the regulation of ke eping cell size in their models ([Bottino et al 



2002; iRubinstein et all l2005| ). Here, we also impose the rule of keeping cell 



size instead of considering the detailed mechanisms of contraction and cortical 
tension. This simplification might be partly justified by the experimental data 
that the shapes of wildtype Dityostelium discoideum does not differ much from 
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those of myosin- II-mill mutants even though the exerted mechanical forces are 
clearly different ( Uchida et al. . 2003f ). The amiB and myosin //double nockout 
mutants are also known to show a ke ratocyte-like shape as in mutants which 



only lacks AmiB ([Asano et all I2004T ) . 



(4) Cellular domain extension: The rule randomly selects a grid from the 
membrane. When Fj at the selected jth grid in the membrane exceeds the 
threshold F^, an external grid in the six nearest grids of the jth grid is turned 
into a cellular grid. F^ roughly defines the amounts of actin filaments per a 
grid. When there are two or more than two external grids around the jth grid, 
a grid is randomly selected. If this grid is referred to as I, Fi = Fj/2 and other 
variables are set to zero. Fj equals to Fj/2 by definition, where the prime 
indicates the value at the next time step. 

(5) Sampling: The event of sampling does not alter the system but the cellular 
shape, position, and concentrations of molecules are monitored. "One step" of 
the cellular dynamics is counted when Sampling is called once. 

We also give a master rule that randomly selects one of the above rules. Prob- 
abilities of selection for rules from (1) to (5) are written as Pi, P2, P3,Pi and 
P5, respectively with Pi + P 2 + P3 + P4 + P5 = 1- When the master rule se- 
lects one of the five rules, the selected rule is executed. We iterate this process 
several millions times. 



We assume that the length of a grid is approximately 1 /^m, we put the initial 
shape of the cell to be a circle with 30-grid diameter, and the equilibrium 
volume is set to be Vq = 900. The effective diffusion constant of the inhibitor 
is P e g = (P/6)(1.0yum) 2 /5t x (P 2 /P 5 )/V r , where 5t is the time length of one 
step. By setting 5t = 0.3s, P 2 = 0.8718, P 5 = 0.00029, and D = 0.45, we 
have P e ff ~ 0.8/xm 2 /,s, which is of the s ame order of the diffusion constant of 



membrane proteins ( Gamba et al. , 20051 ) . Other parameters are set to prevent 



the actin filament from spreading too broadly along the membrane but to 
be hete rogeneously distrib uted in response to the anisotropic environmental 

~T); /3 = 0.1, k a = 0.9, k p = 0.02, 7 = 4.0, k f = 0.99, 
0.040686, P, = 0.05812, P 4 = 0.02906, c = 1.2 and 



stimuli flii ima et al 
h = 10.0, F th = 1.0, Pi 
kT = 20. 



The cell locomotion is simulated with two different patterns of the external 
signal distributions Sj\ One is "linear gradient" Sj = 0.0222?/, + 1, where jjj 
indicates the position of the jth grid along the y direction in the grid space, 
and the other is the uniform distribution Sj = 1.0. The cell locomotion follows 
the gradient of the signal distribution in the former case. Even in the uniform 
signal distribution of the latter case spontaneous polariza tion and unidirec - 
tional movement of cell have been experimentally observed ((Wilkinsonl . 11998) . 
In the present model intracellular fluctuations of activator and inhibitor are 
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amplified by the fluctuation of the cell shape and grow into the spontaneous 
cellular locomotion, for which we make observation on the relation between 
motility and morphology. 



3 Results 



When the parameter a, the rate of activator creation, is set to 1, the long axis 
of the cell shape roughly coincides with the moving direction of the cell (Fig. 
1(a)). The external signal is distributed to have the linear gradient. Fig. 1(b), 
(c) and (d) indicate localization of activator, inhibitor and actin filaments, 
respectively. 

Activator accumulating at the membrane is more localized around the front of 
cell than the rear, whereas inhibitor distribution is biased towards the rear. As 
a consequence of such localization of activator and inhibitor, actin filaments 
are created at the front of the cell. These localized patterns are similar to 
those observed in Dictyostelium discodeum f liiima et all I2OO2I ). where PI3K 



and F- actin are biased towards the front, whereas PTEN accumulates at the 
rear of Dictyostelium, showing that the cell with a = 1 is the neutrophil type. 

When the parameter a is set to 1.6, the shape becomes keratocyte type, that 
is, the short axis of the cell shape roughly parallels its moving direction (Fig 
. 2(a)). The external signal in Fig. 2 has the same linear gradient as in Fig. 
1. Localization of activator and inhibitor again leads to creation of actin fila- 
ments at the front of cell as shown in Fig. 2 (b), (c) an d (d). These features 



agree with the experime ntal observations of keratocytes ([Svitkina et all 1 1997 
Verkhovskv et all 19991 ) and those of the keratocvte-like cells ()Asano et al " 



20041 ). Note that the cellular shape with a = 1.0 or a = 1.6 in the uniform 



external signal is almost same as that in Fig. 2. 

In order to analyze the cell shape quantitatively, we calculate the direction of 
"minor axis" by fitting an ellipse to the the cell shape; 



N 



[&x i &y ) 



' N 
i=l 



M 



(6) 



where e x and e y are x and y-components of a vector pointing along the minor 
axis and Z is the normalization factor to make e x + e y = 1. (xio, yio) is relative 
position of the ith cellular grid measured from the center of mass of the cell. M 



is calculated as M 



N 



1 x i0t m yy 



2 \J^xx ~t~ ^^yy 
£*Li ViO and m xy = 



+ y {TTl xx 
: l^i=l x iOViO- 



m 



mi 



l 2 + Am? xy 



where m 7 



The cellular velocity (v x (t n ),v y (t n )) at time t n is calculated as 
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v x {tn) = (9x(tn+m)-9x(t n ))/mSt and v y (t n ) = (g y (t n+m )-g y (t n ))/m5t, where 
{9x{t n ), 9y{t n )) indicates center of mass of the cellular domain at step t n , and 
5t is a step length in the simulation. We select an integer m so that v(t) does 
not pick up the small fluctuation. 

Correlation between the axis and the cellular velocity, C r , is calculated by 
C r =< \v x e x + v y e y \ > /\f(vl + Vy) >, where the bracket <> indicates average 
over both the time interval in each simulation run and runs started with 
different random seeds. If there is no correlation between velocity and the 
direction of minor axis, C r should equal to ^ J 2 * \cos9\d9 = 2/tt ~ 0.6366. 
Therefore, we define that the cell is keratocyte and neutrophil type if C r is 
significantly larger and smaller than 0.6366, respectively. Fig. 3(a) indicates 
C r as a function of the parameter a, which has a minimum around a — 1.0 
and a maximum around a = 1.6, showing that the cell becomes neutrophil 
type around a = 1.0, and becomes keratocyte type around a = 1.6. Note that 
C r seems to be closer to 0.6366 when a is smaller, showing that the correlation 
vanishes between the cellular shape and the velocity direction. This is because 
the cellular shape becomes approximately spherical and the cell fluctuates 
isotropically as a approaches 0. 

Fig. 3(b) shows the relationship between the averaged cellular speed 



v =< ^Jv^ + Vy > and a. The peak exists around a = 1.6, located at the 
same a at which the correlation C r shows the peak, which implies that the 
keratocyte type cell runs faster than cells with other value of a. 

Fast movement, however, should require more amount of chemical energy. 
Since ATP is consumed in actin polymerization, amount of ATP molecules 
used or amount of chemical energy required should be proportional to the 
number of grids, N g , created by the fourth rule. Fig. 3(c) shows the relationship 
between a and W =< 5N g /5t >, where 8N g is the number of grids created by 
the fourth rule per one simulation step. W gradually increases as a increases. 
If energy source would be limited and a cell dose not need to rush to its 
destination, it is efficient for the cell to reduce "energy per unit length" . We 
estimate "energy per unit length" as follows: 



where T indicates the time interval of each simulation. Fig. 3(d) indicates that 
there is a minimum of E ad around a = 1.0, implying that the neutrophil type 
cell moves most efficiently 

In order to arrive at cell's destination rapidly, the cell should move as straight 
as possible. To estimate this straightness we calculate S t by the distance from 
the start point to the end point divided by the total path length of the tra- 





(7) 
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jectory of the center of mass of the cellular domain. St ranges from to 1 by 
definition. St approaches 1 when the trajectory becomes completely straight. 
Fig. 3(e) shows that in the uniform distribution of the external signal (indi- 
cated by a solid line) a peak exists at a = 1.6, whereas in the linear gradient 
distribution (dashed line) a flat peak exists from a = 1.2 to 1.6. Although there 
is a slight difference in St of keratocyte type (a = 1.6) and St of neutrophil 
type (a = 1.0), straightness of two types seems almost equivalent. 



4 Discussion 





a = 0.4, 7 = 100 


a = 1.0, 7 = 4.0 


a = 1.6, 7 = 4.0 




0.51 ±0.15 


0.53 ± 0.0065 


0.86 ± 0.0030 


V 


0.11 ±0.0028 


0.088 ±0.001 


0.12 ±0.00057 


w 


11 ±0.26 


6.2 ± 0.068 


13 ±0.041 


E a d 


97 ±1.9 


71 ±0.39 


110 ±0.51 


s t 


0.48 ± 0.035 


0.61 ± 0.0085 


0.68 ± 0.0064 



Table 1 



Our model shows that the keratocyte type cell moves more rapidly and more 
straight than other types of cells. These feat ures agree with experimental 
obser yations that kerato cytes ( Lee et al. . 19931 ) or the keratocyte-like amiB~ 



cells (jAsano et all 120041 ) crawl rapidly along a straight direction. In our model 



this rapid motion is because of the wide spreading cell membrane at the cell 
front: Since actin polymerization takes place in the wide area in the keratocyte 
type cell, the driving force of the forward protrusion is stronger than in the 
cell with a narrow leading edge. 

Another way to realize the rapid motion is to accelerate the polymerization 
itself. In our model the actin polymerization rate, 7, controls the rate of actin 
polymerization. We increase the value of 7 to 100, meaning that the actin 
polymerization rate increases twenty-five times. This large change of 7 is con- 
sistent with the fact that the rate of actin polymerization can sensitively vary 
tens or a hundred times according to the co ncentrations of G-actin, profin, and 
other related proteins ( Kovar et al. . 20,06). In order to prevent the cell front 
from expanding explosively, The creation rate of activator, a, is suppressed 
to be 0.4. Table 1 compares C r , v, W, E ad and S t for three parameter sets 
with different 7 and a in the linear gradient condition. The parameter set 
(a = 0.4, 7 = 100) allows the cell to move as fast as the keratocyte type but 
the motion is less straight. The cell with this parameter set is the neutrophil 
type having C r value less than 0.6366; C r = 0.51 ± 0.15 as summarized in 
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Table 1. Its detailed shape, however, is different from that in Fig. 1(a). Large 
fluctuation of the cell shape as shown in Fig. 4 is characteristic to this rapidly 
moving neutrophil type cell, which leads to the less straightforward motion 
of the cell. Consumed energy is less than the keratocyte type. The wild-type 
Dictyostelium in the starved phase moves rapidly but along the fluctuating 
direction and seems to correspond to this type of cell. 

Efficient energy usage of the rapidly moving neutrophil type predicted in the 
model is consistent with the observed decrease in activated NADH in Dic- 
tyostelium discodeum cells: Amount of activated NADH should reflect the 
rate of energy consumption, which is large in unpolrized cells in the vegeta- 
tive phage but decras es in highly polarized elongated cells in the starved phase 
((Nomura et allbOOfih . 



Why does a cell with large 7 and small a have a fluctuated "Dicyosteium"-like 
shape? We can guess the mechanism by following the simulated behavior of 
cell: Small a suppresses actin polymerization, but once the actin polymeriza- 
tion is initiated as a fluctuation, the large 7 leads to a rapid protrusion. In this 
way the large 7 and small a makes protrusion abrupt and intermittent. This 
feature of movement of Dicyosteium-like cell gives a sharp contrast to that of 
the keratocyte type cell. The keratocyte type cell with a larger a but with 
smaller 7 shows the continuous actin polymerization, leading to the smooth 
progression of cell. Thus, the cellular shape and movement decisively depend 
on 7 and a. Here, we can summarize how the cell shape depends on these two 
parameters: Small 7 and small a lead to the slowly moving neutrophil type, 
small 7 and large a lead to the keratocyte type, large 7 and large a lead to 
the explosive expansion, and large 7 and small a lead to a Dicyosteium-like 
shape. It should be interesting to further analyze the morphological change in 
the two-parameter space of 7 and a in a quantitative way. 

We omitted describing detailed molecular processes of contraction at the rear. 
Instead, we assumed the rule of keeping cellular size. Although our results 
showed that this assumption explains several patterns of amoebae, the assump- 
tion clearly does not explain the cell locomotion which shows the temporal sep- 



arati on between the protrusion phase and the contraction phase ([Uchida et al 



2003b . We should improve the model to explicitly take account of the contrac- 



tion processes, which is left as a further avenue of research. 



5 Conclusion 



By using a theoretical model we showed that the rate of production of activa- 
tor is a key factor to determine whether the cell shape becomes the neutrophil 
type or the keratocyte type. The apparent morphological conversion between 
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two types can be caused by rather simple biochemical or genetic perturbations 
which modulate the activator production rate. The keratocyte type cells crawl 
rapidly along the straight direction by consuming a large amount of energy. 
Locomotion of the neutrophil type cell is most economic in terms of energy 
consumption. This slower neutrophil type cell can crawl as straight as the ker- 
atocyte type cell. The neutrophil type cell can crawl as rapid as the keratocyte 
type cell when the rate of actin polymerization is large enough. The moving 
direction of this rapidly crwaling neutrophil type cell is more fluctuating than 
the keratocyte type cell. A plausible explanation of the reason why starved 
Dictyostelium cells take strategy of the rapid neutrophil type locomotion is 
that they have to gather rapidly but the fine tuning of moving direction with 
fluctuating trials and errors should be needed when they approach the aggre- 
gation center. Further systems biological investigations through simulations 
with the model which treats the morphological change and chemical reactions 
in a unified way should enable to classify a variety of different strategies of 
locomotion and should provide a guide line to design experiments. 
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Fig. 1. (a) A snapshot of the cell shape with a = 1.0. An arrow indicates the 
direction of cellular motion (b) Distribution of activator in the same cell as in 
(a). Local densities are expressed with gray scale, (c) Distribution of inhibitor, (d) 
Distribution of actin filaments. 
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Fig. 3. Relationships between the parameter a and five values, C r , v, W, E a< x and 
St. Solid and dashed lines indicate simulations in the environment of the external 
signal having the linear gradient and the uniform concentration, respectively. The 
horizontal axis of all graphs is the parameter a. (a) The vertical axis is C r : cor- 
relation between cellular velocity and the semimajor axis of the cell. Error bars 
indicates standard deviation. All points are values averaged over 50 different sim- 
ulations started from different random seeds, (b) v: Averaged speed, (c) W: The 
number of grids created by the fourth rule per unit time, which is thought to be 
proportional to the number of consumed ATPs. (d) E a d'. the number of grids created 
by the fourth rule per unit length, (e) Sf. The distance from start point to end point 
divided by the total path length of the trajectory of the center of mass. 
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Fig. 4. A snapshot of the cell with a = 0.4 and 7 = 100. 
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